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Abstract  -  This  study  demonstrates  a  technique  of  estimating  soil 
moisture  using  the  combined  passive/active  L-band  microwave 
measurements.  It  shows  1)  evaluation  of  the  small  albedo 
assumption  for  using  dual  polarization  passive  measurements,  2) 
development  of  a  synthesized  technique  to  estimate  soil  moisture, 
and  3)  evaluation  with  ground  soil  moisture  measurements  from 
the  SMEX02  experiment  data. 

I.  INTRODUCTION 

Soil  moisture  is  a  key  parameter  in  numerous 
environmental  studies,  including  hydrology,  meteorology,  and 

agriculture.  It  plays  an  important  role  in  the  interactions 
between  the  land  surface  and  the  atmosphere,  as  well  as  the 
partitioning  of  precipitation  into  runoff  and  ground  water 
storage.  Therefore,  the  spatial  and  temporal  dynamics  of  soil 
moisture  are  important  parameters  for  various  processes  in  the 
soil-vegetation-atmosphere-interface.  The  Hydrosphere  State 
Mission  (Hydros)  with  both  Active/Passive  L-band 
instruments  has  been  approved  by  NASA  for  monitoring 
global  soil  moisture  and  freeze/thraw.  The  Hydros  instrument 
combines  radar  and  radiometer  subsystems.  The  radar 
operates  with  VV,  HH,  and  HV  transmit-receive  polarizations, 
and  uses  separate  transmit  frequencies  for  the  H  (1.26  GHz) 
and  V  (1.29  GHz)  polarizations.  The  radiometer  operates  with 
V,  H  and  U  (third  Stokes  parameter)  polarizations  at  1.41 
GHz. 


In  attempt  to  use  the  active  or  passive  microwave  remote 
sensors  for  estimation  of  soil  moisture,  we  are  mainly  facing 
two  common  problems:  effects  of  surface  roughness  and 
vegetation  cover.  Natural  variability  and  the  complexity  of  the 
vegetation  canopy  and  surface  roughness  significantly  affect 
the  sensitivity  of  backscattering  and  brightness  temperature  to 
soil  moisture.  Backscattering  and  brightness  temperature 
signals  from  the  vegetated  areas  is  a  function  of  water  content 
and  its  spatial  distribution  as  determined  by  vegetation 
structure  and  underlying  surface  conditions  including  surface 


roughness  parameters  and  dielectric  properties.  Due  to  the 
limited  observations  from  either  passive  or  active 
measurements  alone,  an  ill  condition,  the  number  of 
measurements  and  equations  are  less  than  the  number  of 
unknowns,  is  expected.  It  results  in  the  uncertainties  in 
estimation  of  soil  moisture. 

In  this  study,  we  study  a  combined  active/passive 
technique  to  estimate  surface  soil  moisture  with  the  focus  on 
the  short  vegetated  surfaces.  We  first  simulated  a  database  for 
both  active  and  passive  signals  under  Hydros’s  sensor 
configurations  using  the  radiative  transfer  model  with  a  wide 
range  of  conditions  for  surface  soil  moisture,  roughness  and 
vegetation  properties  that  we  considered  as  the  random 
orientated  cylinders.  Using  this  database,  we  developed  1)  the 
techniques  to  estimate  surface  soil  moisture.  We  will 
demonstrate  this  technique  with  the  model  simulated  data  and 
its  validation  with  the  airborne  PALS  image  data  from  the  soil 
moisture  SMEX’02  experimental  data. 

II.  INVERSION  MODEL  DESCRIPTION 


The  simulation  models  used  in  this  study  are  based  on 
the  commonly  used  three  components  emission  and 
backscattering  model  concept  with  no  consideration  on  the 
polarization  dependence  of  the  vegetation  extinction 
properties.  They  are  given  as 


Tbp  —  Tv  •  (Ep  +  EpV )  +  Ts  •  E*  •  L 


t  _  _v  ,  sv  ,  s  r  2 

°  pq  _  °  pq  “l"  °  pq  0  pq  '  ^pq 


(i) 


The  superscrip  t,  v,  s,  and  sv  indicate  the  total  and  components 
from  vegetation,  soil,  and  soil-vegetation  interaction  terms, 
respectively.  T  is  the  temperature  of  vegetation  and  Ts  is  the 
soil  temperature.  Ep=(l-co)(l-Lp)  is  the  emissivity  of  the 
vegetation  layer.  Epv=  EpRpLp  is  the  emissivity  due  to 
vegetation-surface  interaction.  Lp=exp(-Tp/cos9)  is  the  one  way 
attenuation  factor  or  transmittivity  of  vegetation  layer.  The  x  is 
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the  optical  thickness.  The  direct  volume  backscattering 
component  is  given  by 

O'pq  (#)  =  0.5icfq  (l  -  Lpq )/  xf  q  (2) 

The  surface  emission  and  backscattering  are  simulated  by  the 
Advanced  Integral  Equation  Model  [1-2]. 

III.  EVALUATION  OF  SMALL  ALBEDO  ASSUMPTI  ON  IN 
PASSIVE  MEASUREMENTS 

For  convenience  of  analyses,  we  re-write  (1)  for  passive 
emissivity  as 

E|,  =  /-<y-(/-Lp)-Rp-Lp-[/-EpJ 

=  i-Vc-R*-Ve  (3) 

Where  Vc  =  co  ■  (7  -  Lp  )  is  the  vegetation  correction  factor. 

Ve=Lp-(f-Ep)  *s  vegetation  attenuation  factor.  Due  to 

Vc  typically  quite  small,  the  study  showed  that  the  effect  of 
Vc  wss  small  on  estimation  of  soil  moisture  when  using  H 
polarization  and  was  assumed  that  co=0  [3].  This  is  because  1) 
the  vegetation  effects  in  low  frequency  passive  signals  are 
mainly  dominant  by  the  absorption  and  2)  the  surface  effective 
reflectivity  in  H  polarization  has  much  larger  impact  on  the 
measured  emission  signals  than  that  from  the  scattering 
generated  signals.  With  the  small  albedo  assumption  of  co=0, 
(3)  becomes 

Ep  =  7  —  Rp  ■  Lp  (4) 

To  evaluate  this  approach  on  estimating  soil  moisture  for 
using  the  dual  polarization  measurements,  we  simulated  the 
emission  signals  using  (3)  with  the  t:  0.01  to  0.6,  co:  0.01  to 
0.05  and  a  wide  range  surface  soil  moisture  and  roughness 
conditions.  Figure  1  (left)  shows  the  histogram  of  Vc  from 
this  simulation.  It  can  be  seen  that  Vc  are  quite  small  with  a 
range  less  than  0.03.  In  order  to  examine  its  impact  on  the 
description  of  the  relationship  between  V  and  H  polarizations, 
we  compared 

^  ~  Vc  ~  Ey  _  R^Ve  _R°  (5) 

1  -  Vc  -  Ej,  R|  •  Ve  Rh 

Figure  1  (right)  shows  the  left  side  of  (5)  as  x-axis  that 
represents  the  case  without  getting  co=0  and  the  right  side  of 
(5)  as  y-axis  that  represents  the  effective  reflectivity  ratio 
value  with  assumption  of  co=0.  It  shows  that  the  small  albedo 
assumption  will  result  in  the  significant  effects  on  deriving  the 
ratio  of  the  surface  polarization  reflectivity,  which  will  greatly 
impact  on  the  estimation  of  soil  moisture.  This  is  because  the 
emissivity  of  V  polarization  can  be  quite  close  to  unity  so  that 
Vc  can  be  at  the  same  magnitude  as  the  effective  reflectivity. 


Therefore,  the  effect  of  Vc  has  to  be  take  into  account  in 
estimation  of  soil  moisture  when  using  the  ratio  approach. 
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Fig.  1.  The  histogram  of  the  simulated  Vc  at  left  and  the  comparison  on  the 
polarization  ratio  signals  with  (right  side  of  (4))  and  without  (left  side  of  (4)) 
the  small  albedo  assumption. 

IV.  ESTIMATION  OF  THE  VEGETATION  CORRECTION 
FACTOR  Vc  WITH  ACTIVE  RADAR  MEASUREMENTS 

Natural  variability  and  the  complexity  of  the  vegetation 
canopy  and  surface  roughness  significantly  affect  the 
sensitivity  of  radar  backscattering  to  soil  moisture. 
Backscattering  signals  from  vegetated  areas  is  a  function  of 
water  content  and  its  spatial  distribution  as  determined  by 
vegetation  structure  and  underlying  surface  conditions.  It  has 
been  realized  that  the  different  radar  polarization 
measurements  has  the  different  sensitivities  to  the  different 
surface  properties.  Especially,  the  radar  cross-polarization 
measurements  are  very  sensitive  to  the  vegetation  information 
since  the  surface  backscattering  does  not  generate  significant 
cross-polarization  signal.  As  shown  in  [4],  the  direct  volume 
backscattering  component  of  the  co-polarizations  in  (1)  for  the 
random  orientated  short  cylinders  can  be  directly  estimated 
from  the  cross-polarization  signals.  That  is 

cr^p  =  3  -  crfh  (6) 

In  comparison  of  Vc  and  <7^  ,  it  can  be  found  that  they  are 

highly  correlated.  Their  relationship  can  be  well  characterized 
as 

log(Vc )  =  a  +  b  •  log(Cpp )  (7) 
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Fig.  2.  Comparison  of  of  Vc  of  the  simulated  by  (3)  (x-axis)  and  the 
estimated  from  the  direct  volume  backscattering  component  (y-axis). 

Figure  2  shows  the  comparison  of  Vc  of  the  simulated  data 


using  (3)  and  the  estimated  from  the  direct  volume 
backscattering  component  of  the  simulated  radar  signals  using 
(7).  It  demonstrated  that  Vc  could  be  estimated  from  radar 
cross-polarized  backscattering  measurements  by  using  (6)  and 
(7).  It  improves  the  accuracy  in  describing  the  relationship 
between  V  and  H  polarization  signals  in  the  soil  moisture 
retrieval  by  the  passive  measurements. 

V.  A  SYNTHESIZED  INVERSION  TECHNIQUE  TO  ESTIMATE 
SOIL  MOISTURE  USING  PASSIVE  MEASUREMENTS 

With  the  estimated  Vc  parameter  from  radar,  we  can  now 
obtain  the  measurements 

^-Vc-Ep  =  Rp  Ve  (8) 

For  estimation  of  soil  moisture,  we  derived  a  synthesized 
technique  that  consists  the  three  formulas  with  the  different 
sensitivity  to  the  different  surface  properties. 

The  first  formula  is  the  polarization  ratio  to  derive  the 
effective  reflectivity  ratio  measurements. 

Ry  •  Ve  _  Ry  _  rv  '  Hy  (9) 

Rh-Ve  Rch  rhHh 

where  Hp  represents  the  effect  of  the  surface  roughness  at 
polarization  p.  Since  it  is  assumed  that  there  is  no  polarization 
dependence  for  the  vegetation  extinction  properties  in  this 
study,  the  ratio  measurement  will  cancel  out  the  vegetation 
attenuation  effect  -  Ve.  This  measurement  provides  the 
information  on  the  soil  moisture  and  roughness.  It  is 
insensitive  to  vegetation  properties. 

The  second  formula  is  obtained  from  our  bare  surface 
soil  moisture  retrieval  study  [4].  For  bare  surface,  the  Fresnel 
reflectivity  ratio  can  be  estimated  by 

—  =  exp[a+bTog(Ry)+cTog(Rh)  +  dRy /R|] 

rh 

The  parameters  a,  b  and  c  are  determined  through  regression 
analyses  from  our  A1EM  model  simulated  database.  These 
parameters  are  dependent  only  on  incidence  angle.  This 
model  uses  the  different  weight  on  the  surface  effective 
reflectivity  measurements  at  different  polarizations  to 
minimize  the  surface  roughness  effects  so  that  the  response  of 
the  measurements  to  the  Fresnel  reflectivity  ratio  (left  side  of 
(10))  can  be  estimated.  Under  present  of  the  vegetation,  it 
becomes 

^■Veb+C  =  exip[a  +  b  log(R^ -Ve)+c  log(Reh-Ve)+d  / Reh]  (11) 

rh 


between  the  Fresnel  reflectivity  at  V  and  H  polarizations.  It 
can  be  well  described  as 


r 


v 


which  leads  to 


(12) 


ky  '  Ve  _  Hy  y/~a  (13) 

(Rh  ■  ve)a  _  H| 

from  (8)  .  This  measurement  provides  the  information  on  the 
surface  roughness  and  the  vegetation  attenuation.  It  is 
insensitive  to  soil  dielectric  properties. 

However,  the  roughness  presentation  in  (13)  does  not 
directly  related  to  the  roughness  presentation  in  (9).  In  order  to 
estimate  soil  moisture,  we  need  to  characterize  the  relationship 
between  the  roughness  parameters  in  (9)  and  (13).  Figure  3 
(left)  shows  the  relationship  between  Hv/Hh  as  x-axis  and 
Hv/Hha  as  y-axis  in  our  A1EM  model  simulated  database.  It 
shows  that  these  two  roughness  properties  have  a  fairly  good 
relation.  By  adding  the  effective  reflectivity  ratio 
measurements,  their  relationship  can  be  well  described  as 


fj  fj  p  e 

^  =  f  +  g*_v  +h*_^ 

Hfl  Hh  Re 


(14) 


where  f  g,  and  h  are  coefficients  determined  from  the 
regression  analyses.  Figure  3  (right)  shows  the  comparison 
between  the  left  side  and  right  side  of  (14). 


Fig.  3.  The  relationship  between  Hv/Hh  as  x-axis  and  Hv/Hh“  as  y-axis  (left) 
and  the  comparison  between  the  left  side  and  right  side  of  (14). 

With  above  formulas  (8)  to  (14),  we  are  now  able  to 
estimate  the  Fresnel  reflectivity  ratio,  the  surface  roughness 
ratio,  and  the  vegetation  attenuation  factor.  We  first  use  the 
radar  cross-polarization  measurement  with  (6)  and  (7)  to 
estimate  the  vegetation  correction  factor  Vc.  Then,  we  loop 
through  the  possible  vegetation  attenuation  factor  Ve  and 
select  the  minimum  error  from  examining  (9),  (11),  (13),  and 
(14)  to  find  the  estimations.  After  we  obtain  the  estimation  of 
the  Fresnel  reflectivity  ratio,  it  can  be  converted  to  soil 
moisture. 


This  measurement  provides  the  information  on  the  soil 
moisture  and  the  vegetation  attenuation.  It  is  insensitive  to 
surface  roughness  properties. 

The  third  formula  is  obtained  based  on  the  relationship 


Figure  4  shows  the  histogram  of  the  absolute  errors  on 
estimation  of  soil  moisture  from  the  performance  of  the  above 
technique  using  our  simulated  emission  database.  This 


database  covers  a  wide  range  of  vegetation  extinction 
properties  and  surface  dielectric  and  roughness  properties  with 
total  240,000  simulated  data.  The  RMSE  of  the  estimated  soil 
moisture  is  0.03  or  3  %  for  volumetric  soil  moisture. 

lj 


-0.2  0.0  0.1  0.2  0.3 


Fig.  4.  The  histogram  of  the  absolute  soil  moisture  errors. 

VI.  EVALUATION  WITH  PALS  MEASUREMENTS 

PALS  is  a  non-scanning  real  aperture  combined 
microwave  radiometer  and  radar,  operating  at  1.41  and  2.69 
GHz  (radiometer)  and  1.26  and  3.15  GHz  (radar)  with 
multiple  polarizations.  The  instrument  was  designed  for  high 
accuracy  measurements  of  ocean  salinity  and  soil  moisture. 
The  radiometer  operates  at  V  (vertical)  and  H  (horizontal) 
polarizations,  while  the  radar  operates  at  VV,  HH,  and  VH 
polarizations.  For  SMEX02  this  angle  was  fixed  at  45°.  The 
instrument  thus  samples  a  single  footprint  track  along  the 
flight  path.  At  a  nominal  flight  altitude  of  1  km  and  incidence 
angle  of  45°  the  instantaneous  3-dB  footprints  at  the  surface 
are  approximately  330  x  470  m.  Additional  details  of  the 
instrument  design  and  engineering  data  are  provided  in  [4]. 

PALS  flew  several  flight  lines  each  day  over  the  test 
area  from  June  25  to  July  8,  2002.  During  the  late  June,  the 
soil  was  very  dry.  At  early  July,  there  were  several  rainstorms 
moved  through  the  area  providing  an  opportunity  to  observe  a 
large  dynamic  soil  moisture  range.  Vegetation  types  in  the  test 
area  are  mainly  corn  and  soybean  with  the  field-averaged 
vegetation  water  contents  observed  on  the  ground  primarily 
within  the  a  few  kg  m'2  range  during  the  SMEX02  experiment. 
The  ground  soil  moisture  at  different  soil  depth  was  sampled 
over  32  field  test  sites  intensively  during  the  experiment. 
There  are  total  of  182  soil  moisture  data  available  with  the 
corresponding  PALS  measurements  from  6  different  days. 

Figure  5  (left)  shows  the  comparison  between  the 
averaged  volumetric  soil  moisture  from  the  ground 
measurements  and  that  estimated  using  L-band  PALS 
measurements  by  the  technique  described  in  last  section.  In 
overall,  the  estimated  soil  moisture  matches  the  trend  of  that 
measured  on  the  ground  and  with  a  reasonable  accuracy.  The 
accuracy  in  terms  of  RMSE  is  0.047  or  4.7  %  of  the 
volumetric  soil  moisture.  Figure  5  (right)  shows  the  histogram 
of  the  absolute  errors  of  the  estimated  volumetric  soil 


moisture.  As  we  can  see  that  there  is  a  system  error  or  bias  in 
the  estimation.  The  largest  error  can  reach  10  %.  Through  our 
analyses,  it  was  found  that  the  estimation  for  the  soybean 
fields  had  a  better  accuracy  than  that  from  the  corn  fields. 
This  may  be  resulted  from  our  assumption  on  the  no 
polarization  dependence  of  vegetation  cover  in  our  inversion 
technique.  However,  corn  has  a  strong  preferred  special 
structure  and  orientation.  It  results  in  a  significant  polarization 
dependence  in  the  extinction  properties  of  the  vegetation, 
which  has  a  great  impact  on  soil  moisture  estimation. 

Fig.  5.  (left)  Comparison  of  the  estimated  soil  moisture  from  the  ground  in- 
si  tu  measurements  (x-axis)  and  that  estimated  from  the  PALS  L-band 
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measurements  (y-axis).  (right)  The  histogram  of  the  absolute  errors  of  the 
estimated  volumetric  soil  moisture. 

VII.  Summary 

In  this  study,  we  develop  a  combined  active/passive 
technique  to  estimate  surface  soil  moisture  with  the  focus  on 
the  short  vegetated  surfaces.  It  shows  that  the  small  albedo 
assumption  will  result  in  the  significant  effects  on  deriving  the 
ratio  of  the  surface  polarization  reflectivity,  which  will  greatly 
impact  on  the  estimation  of  soil  moisture.  We  developed  a 
technique  to  estimate  soil  moisture,  which  includes  the 
estimation  of  the  vegetation  parameter  of  Vc  using  the  cross¬ 
polarization  measurements  of  radar,  a  synthesized 
measurements  from  radiometer.  The  test  from  the  SMEX02 
PALS  measurements  performed  reasonable  well  with  the 
RMSE  of  4.7  %  for  volumetric  soil  moisture  estimation.  It  is 
also  found  the  polarization  dependence  of  the  vegetation  cover 
has  a  great  impact  on  the  soil  moisture  estimation  and  needs 
to  be  considered  in  future  study. 
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